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QCD radiative correction to pair-annihilation of spin-1 bosonic Dark Matter 
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The next-to-leading order (NLO) QCD corrections are calculated for the pair-annihilation of 
spin-1 dark matter (DM) by dimensionally regularizing both ultraviolet and infrared singularities 
in non-relativistic limit {v <^ 1). The complete 0{as) correction is about 8% due to the massless 
gluon contribution. An extra 5% will be added if there is a new interaction from a massive gluon 
of approximately same mass as the DM particle. The NLO QCD correction could give the sizable 
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,.^ ' shift to the DM mass constrained by relic density measurements. 
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I. INTRODUCTION 



I ^ The spin-1 dark matter candidate B' , which is the Z2-odd partner of the hypercharge gauge boson B, appears in 

> 

^\ ' interesting models like universal extra dimensions (UED) [1] and Little Higgs (LH) [2], motivated by solving the gauge 

■^ ' hierarchy problem. The phenomenology of spin-1 DM (B') has been well studied at the tree level. With advances 
f — in precision measurements in cosmology and astrophysics as well as the advent of LHC era, the NLO calculation 

o 

g ot the B-B' annihilation bccon>.s tin.ely f„ accu,-... .n.ly.si.,. We p.™. .he „ex.-.„-lea<ii„g order (NLO) QCD 
^ , corrections for the pair- annihilation of spin-1 bosonic dark matter by dimensionally regularizing both ultraviolet 
(UV) and infrared (IR) singularities in the non-relativistic limit {v ^ 1). Our analysis of the process {B'B' -^ qq) 
is applicable to general spin-1 DM, as in UED, LH, etc. We assume that the quark field q interacts with its Z2-odd 
partner q in the form. 



C D --gyqi'^iYLPL + YRPR)qB'^ + h.c. (1) 

The co upling gy and the number Yl in the UED would be the usual hypercharge coupling and the corresponding 
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FIG. 1: The Born diagrams for B'B' -^ gg. The bold wavy hnes represent the DM candidate gauge bosons {B'), the bold 
solid lines represent the heavy quarks (the Z2-odd partners of the outgoing quarks) and the light solid lines correspond to the 
quarks. 



quantum number, 



9Y = TZ^^ YL{u,d) = l, Yr{u) = 1, YRid) = -l 



cos 9w 



6' 



However in the Littlest Higgs model (L^H) [2], YL{u,d) = -^ and otherwise Yr ~ 0. 

We also assume that the mass M of g is not too far above the mass M of the spin-1 DM particle B'. The approximate 
relation M ss M is valid in UED, but a choice of parameters in other models is required. However, such a choice 
allows us to carry through analytical calculation and give simple results. 

Fig.l shows the Born diagrams for B'{pi)B'{p[) -^ q{p2)q{P2) ■ The exchange of the Z2-odd quark line is bold-faced. 
The amplitude for the Born diagrams is given by 



Mb = -91-YMP2) [r ^\ ^jj^ 7- + I'^^^-Ai^A PLv{p'2)e,{pi)eM) + RH (2) 



t- IVP 



u- M^ 



In the extremely non-relativistic limit (w ~ 0), pi ~ p'l ~ (p2 + P2)/2 — (-^7 0)- The invariant Mandelstam 
parameters are 



t={-pi+P2? = {p'l-p'2?^-M^ 
U={-p[+P2f = ipi-p'2f^-M^ 



(3a) 
(3b) 
(3c) 



The formula reduces to 



2 ^2^,„ ^ (P2 - p'2Yl'' + (P2 - p'^Yl^ 



Mb - glYMp^ Y^' \' - " PLv{p'2)e,{piYAp'i) + RH (4) 



The consistent annihilation rate in d = 4 — 2e dimensions is 



(1 - 6)r(l - e) /4V Y 2g^(y^4 ^ Y^)N, M^ 



r(2-2e) \U'P ) 



9tt 



(m^ + mA 



where Nc is number of colors, v denotes the relative velocity between spin-1 dark matter candidate pair, ^ is an 
arbitrary mass scale and as is the Born cross section. 

II. NEXT-TO-LEADING ORDER CALCULATION 

The one-loop amplitudes consist of the virtual corrections from the diagrams (a),(b),(c),(d) of Fig. 2. Adding up 
those contributions must be ultraviolet finite. The amplitudes contain infrared divergences due to massless gluon 
virtual exchange. The infrared should cancel exactly against the one present in the gluon final state radiation, 
Fig.3(a),(b),(c). The method is to use dimensional regularization in d = 4 — 2e dimensions with massless on-shell 
quarks to regularize both types of divergences, UV and IR divergences. The individual diagrams are calculated in the 
renormalizable Feynman gauge (the gauge parameter, ^ = 1), which provides the gluon propagator. All new particles 
involved in the calculation are set up to have the same mass as B' . 

A. Virtual corrections 

All the virtual corrections could be expressed by two form factors^ related to tensors, {p2 — p'-i)^^^ + (P2 — p'^fl^ 
and (p2 +^2)^7''' + (-P2 -^p'-i)^!^-, for the massless outgoing particles in static limit. However (p2 —p'-i)^!^^ + (P2 —p'^Yl^ 
only survives, p^ + p'2 is a timelike vector and B' polarization vectors, ^^(pi) and e^(p']^), are spacelike. The tensor 
(%>■? + Vo )'^^'^ + {P2 + P^T^ disappears as contracting the B' polarization vectors. So the virtual corrections are 

^ Two more form factors are possible, which are related to the tensors {^2 i ^'2)9^^ ■ However both do not give any contribution for the 
massless outgoing particles. 
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FIG. 2: The Feynman diagrams which contribute to the next-to-leading order (NLO) QCD virtual radiative correction(_B'B' —f 
qq). The crossed (w— channel) diagrams are not displayed. The bold wavy lines represent the DM candidate gauge bosons {B'), 
the bold solid lines represent the heavy quarks (the Z2-odd partners of the outgoing quarks), curly lines represent the gluons 
and the light solid lines correspond to the quarks. 



expressed with only one form factor, which is the coefficient of the Born amphtude. 



Ml = Fb'Mi 



(6) 



We use the conventional approach of Fcynnian parameters to calculate the corrections. The virtual corrections are 
simplified in the common integral with the shifted momentum and Feynman parameters. The infrared divergences 
appear in the Feynman parameter integrations. 






i{-iy—m Y(r + rf/2)r(m - r - d/2) 



(27r)° 



{P - C)™ (477)^^/2 



r(rf/2)r(m) 



dx,5 {Y.X, - 1) C'-™+'^/2 (7) 



The amplitude of the diagrams in Fig. 2 (a) is 



,(a) 



MTi - -gi^Yiuip,) 



pM , h~ h + M y ^ ,, ^ ^_^ I ^2 



(a) 



t- M^ 



J +T 



M 



t-M"^ 



ria) + ^^^'^ 



PLv{p'2)e^{pi)e^{p'i) (8) 



r^^N and F/'^s are the vertex corrections in rf— dimensions, which are given by 



(a) *^«^^y (27r)'i F(p2+fc)2((p2-pi + fc)2-M2) ^"""^ 



where Qs is the strong couphng and Cp is the Casmir operator of the fundamental representation in the color group. 
The Lorentz indices /i and v are just switched for the ?i— channel, since the invariant Mandelstam parameters t and 
u are identical and the momenta, pi and p'j^, are timelike in the extremely non-relativistic case(i' ~ 0). 
The form factor of Fig. 2 (a) results in 



where the subscript UV implies the 1/e pole coming from UV divergence and Q!s(= gl/ATi) is QCD fine structure 
constant. 

Calculation^ of Fig. 2(b) requires careful and tedious effort because of four propagators in the loop and a three 
folded integral over Feynman parameters. The scattering amplitude is 



M^ix - -9l-Ylu{p2)Tl^^PLv{p',),,{p,)e.{p',) (11) 



Y^l^\ is the combined vertex correction of the t— and m— channel diagrams. It is given by 



pp.^_.2^ ( d^k Yih + /^)7^(i'2- h + h + M)7"(-i>2 + ^)lp + M ^^ i^ . . 

W '^' ^ J {2^Y k^p2 + kfi-p'^ + k)mp2-pi + ky-M^) ^ ' 

This integration could also be manipulated into the common form, Eq.(7), however the integral has singularity 
in the euclidean region. The imaginary parts are included to continue this integral in the euclidean region for the 
positive s and it results in a complex form factor. 

^ The analogous calculation for the box diagram is in Rcf . [3] for the different phenomenology. 



The imaginary parts are not relevant for real radiative nieseurenients as a consequence of the unitarity of the 
S- matrix: 5^5* = 1. 

The contribution of Fig. 2(c) comes from propagator corrections to on-shell quark lines. For the massless quark 
lines, there is no contribution using dimensional regularization since the same regulator is used. In the leading order 
of e, it can be written as 



^(^' - "-^ f ^ V r(i + .) f-^ + ^ 



Fi,''-^^ri&]n^+'>{-^+^\ (») 



However the contribution of Fig. 2(d) comes from off-shell heavy quark lines. It produces a contribution without IR 
divergence. 






F^^>^2--^\,.^_^. = ^^^{--fJ^] r(l + .)(-,^-l + log2) (15) 



where the factor of 2 comes from the two vertices and those give identical contributions. The results show that UV 
divergences are exactly canceled. 

Adding up all the virtual corrections, the QCD corrections are 



.(virutal) „p .p ^ asCp f^TTn^y , nA 1 ^ 27r2 14 8 \ 

B. Real corrections 

The real QCD correction appears in the ratio of annihilation rates for two and three body final states. In an average 
over the polarization of the incoming vector bosons and a sum over the spin and color of outgoing quarks and gluon, 
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FIG. 3: The Feynman diagrams which contribute to the next-to-leading order (NLO) QCD real radiative correction(_B'B' -^ 
gqq). The crossed (m— channel) diagrams are not displayed. The bold wavy lines represent the DM candidate gauge bosons 
(B'), the bold solid lines represent the heavy quarks (the Z2-odd partners of the outgoing quarks), curly lines represent the 
gluons and the light solid lines correspond to the quarks. 



the annihilation rate for three body final states is 



1 



4M2 9 



^Jd^3\Ma+Mb+Mc\^ 



(17) 



Ada, -Mb and Aic are the scattering amplitudes corresponded to the diagrams (a),(b),(c) of Fig. 3, and are given by 



Ma,L = 



. 2 r;2^a"fe)7''(i'2 + /:)(7^i'a7'^+7'^A7^)^L«(p'2)eM(Pi)e.bi)er(fc) 
-i9y9sYlT 



-iglgsYlT 



{p2 + kf{pl-M^) 
a nP 2) irhr+Yhr) (-/2 - /ch''PLv{p'2)ef,{pi)e,ip[)e;'^ik) 
i-p'^-kYipl-M^) 



. 2 .72^a ^(P2) [7^ A7^ A7" + M^-jl^YY + M ^ ^] PL^(p^)£M(pi)£.(rf)£;°(fc) 



(18a) 
(18b) 
(18c) 



where T" is the QCD generator and pa = p'l — P2 = —pi + P2 + k,ph = —pi + p2 = p'l — p'2 — k in shorthand where 
piand p[ are momenta of the incoming vector bosons, p2 and P2 for the outgoing quarks and k for the radiated gluon. 
Pi ~ p'l ~ (p2 + P2 + k)/2 ~ (A'/, 0) in the non-relativistic limit and the numerator of A4c is allowed to express by 
the momenta, p — P2 ^ p'2 and k on the diagrammetic symmetry. 



K - P'ip^l" + P^r) + iPaks^^e^P'^igtY + gir - .9^"7a) + {P2 ■ p'2 + 1M^){g'''l'' + .g"''7^ - g^^^Y) (19) 



The relation, ep(p)e*(p') = —gpa, is adopted to calculate the squared amplitude for both incoming massive spin-1 
dark matter candidates in Feynman gauge, since the polarization vectors efi{pi),e^{p'i) are spacelike. The soft and 



coUinear singularities appear in the squared amplitudes, |A^ap,|A^;,p and 2Re(A^*A^b). 1/e poles are canceled out 
in 2Re(Al*Alc) and 2Re(Al^Alc) as combining the numerator and denominator, and \Mc\ is totally free of IR 
divergences. 

The new dimensionless parameters are introduced for the phase space integration by energy-momentum conserva- 
tion. 






1/e poles appear at xi = and X2 ~ 0. The Lorentz-invariant three body phase space takes the form with the new 
parameters. 



4M2 U-kv?\ 



^'^^^ " (47r)3r(2 - It) \Tm^) Idxid^2dx3{xiX2X3) 'Si^x, ~ 1) 
The overall real correction gives 



(21) 



,(.eal)_«^/4VYr. .M J_ 397 \ 



When all the contributions are added, the infrared are exactly canceled and the finite result is 



SQ^^^—-l- + -log2--] (23) 



The correction is about 8% enhancement for Cp — 4/3 and ^^(ITcV) ~ 0.09, and the finiteness of this correction 
implies t hat there is no divergence by degeneracy. 

^ The corrections can also be approached by the optical theorem to acquire the cross section. But eventually both are identical and we 
would have the same QCD correction. 



C. The corrections by the heavy gluon 

If the heavy Z2-odd of the usual gluon is not far above M, its contribution to the B' pair annihilation could be 
important. The relevant Lagrangian of this heavy gluon field G'^ is 



C D -g^qi^T'^G^q + h.c. (24) 

For being able to carry through analytical calculation, we set the mass of this heavy gluon field GJf equal to the B' 
mass M in our calculation. Such a simplification is valid at least in UED. The corresponding Feynman diagrams are 
similarly given in Fig. 2, except that we need to interchange the bold quark line and the light quark line in the loop 
as well as the usual gluon is replaced by the heavy gluon. Note that there is no IR divergence for a massive gluon, 
and the UV divergence cancels among diagrams. As a single heavy gluon cannot be produced by Z2-symmetry, the 
diagrams in Fig. 3 are all absent and it results in the improved QCD correction. The QCD corrections by the heavy 
gluon can also be calculated fully analytically and the overall correction is 



6V3 



^(gu^on; ^ -- ( _5 + :l^ _ ioi„g2 + ^ ) (25) 



The correction is about 5% enhancement and is comparable to massless gluon correction. 

III. APPLICATION TO RELIC ABUNDANCE 

The spin-1 dark matter candidate B' , which is the Z2-odd partner of the hypercharge gauge boson _B, has been 
an attractive dark matter candidate in UED and LH, but both models were extended in different ways. The spatial 
dimensions are enlarged for UED, otherwise the symmetry groups are enlarged for LH in 4-dimensions. So the masses 
of the new particles are scaled by the extra dimension of size R compactificd on an S^ /Z2 orbifold for UED and 
the enlarged global symmetry breaking scale / for LH. Since the new particles are in the different symmetry group 
structure, the different gauge charges are assigned to the fermions {qqB' coupling) and it causes them to induce the 
different phenomenological analyses. We consider relic abundance with our QCD correction, but the detailed analyses 
are not mentioned since those were well studied in the leading order in the other articles [4, 5]. 
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FIG. 4: Prediction for relic abundance as a function of the WIMP mass. NLO QCD radiative corrections are considered on 
the solid lines. The horizonal band denotes r2h^=0.111±0.012 and defines the WIMP mass window. The WIMP mass in the 
window is shifted about 50 GeV for UED with NLO QCD correction. 



Assuming that the spin-1 B' accounts for all DM relic abundance, one can constrain its mass M based on the pair 
annihilation rate in the early universe. In UED models, the annihilation channel into a fermion pair B' B' -^ ff 
dominates because of sizable couplings of the gauge charges Yj^m. The measured relic abundance prefers M to be 
about TeV. On the contrary, the L^H model has rather small couplings of the gauge charges Yl — -^Xr — 0, and 
the relic abundance only requires M of order of 100 GeV. The detailed quantitative analyses for relic abundance by 
this type of WIMP annihilations can be found in Rcf. [4] for UED and Ref. [5] for LH. Fig. 4 shows the prediction for 
relic abundance as a function of WIMP mass with the present WMAP precision [6] . The WIMP mass in the window 
is shifted about 50 GeV compared to the leading order for UED^, but it shows no difference from the leading order 
by the negligible annihilation fractions into quarks for LH . 

IV. CONCLUSION 



The next-to-leading order (NLO) QCD corrections are calculated for pair-annihilation of spin-1 bosonic dark matter 
by dimensionally regularizing both ultraviolet and infrared singularities in the non-relativistic limit. The order a^ 
correction amounts to about 8% and can enhance to 13%, when including heavy gluon. The NLO QCD correction 
could giv e the sizable shift to the DM mass constrained by relic density measurements. 

** It was truncated to the first level Kaluza Klein (KK) mode for our calculation and UED with such truncation is renormalizable though 
the theory is not renormalizable. 
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APPENDIX: USEFUL FORMULAS AND IDENTITIES 

The Dirac algebra in d = 4 — 2e dimensions are listed in Ref. [7] and they could be induced by the anticommutator 
{7^, j"} — 2g^^^ and the identity 7^*7;^ = d. The regulator e is omitted on trace algebra Tr[I] ~ 2'^/'^ for two and three 
body cross section calculations. 

1/e/ij poles are extracted by the partial integrations on the loop calculations in case of a need to extract the poles 
and the remnants are expanded in ordinary Taylor series with respect to e. The three folded parametric integrals for 
the box diagram are calculated by Feynman parameter properties, since one of them is not symmetric with the others. 



dx dyf{x,y)^ / dy dxf{x,y) (A.l) 



The complicated integrations are avoided and number of integrations is reduced with this property. 
The useful identity to split the singular and non-singular parts for the real corrections is 

TiS^ + r^ (*-^' 

x(L — X) X 1 — X 

This identity can be extended to the higher powers of the variables and will simplify the calculations since the 
regulator e can be dropped for non-singular parts. 

Some of the parametric integrals contain the dilogarithm (or spcncc) function. 



Li,(.) ^ - TrfiMl^i:!) (A.3) 

Jo t 



The values which appear in our calculations are 



12 



Li2(l) = 



6 



Li2(-1) = - 



12 



Lio 



log' 2 
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(A.4a) 
(A.4b) 
(A.4c) 



The parametric integrals which appeared in the box diagrams reduce to the incomplete beta function. 



1/2 



Bi/2(TO,n)=/ dxx"-^ (1 - x) 

"'0 



n-l 



(A.5) 



Bi/2{m, n) + Bii2{n, m) == B{m, n) 



Bii2{m,m) — —B{m,m) 



Bi/2{m,m + 1) = 
Bi/2{m,m + 2) = 



-2m' 



B{m, m + 1) 
B{m,Tn + 2) 



m 

o — 2?n' 



?Tl 



(A.6a) 
(A.6b) 
(A.6c) 

(A.6d) 



where 



B{m, n) 



T{m)T{n) 
T(m + n) 



(A.7) 



For integer n, the Taylor expansion is adopted with respect to e. 
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Bz{-e,0) = 1 + elog{l - z) - e^U2{z) (A.8a) 

B,{-e, 1) = 1 + e (-J^ + log(l - ^)) + e' (log(l - z) - Li2(z)) (A.8b) 



with z = 1/2. 

The parametric integrations involved the heavy gluon loops are 



/ dx\og{x^ -x + 1) = 2(^11-1 

/ dxix"^ - X) l0g(x2 -x+l) = — - -— TT 

Jo 18 6 



\ x^-x 2 

dx^ = —1 



x^ -x + 1 3V3 
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